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Abstract 

Two-phase composites with non-overlapping inclusions randomly 
embedded in matrix are investigated. A straight forward approach is 
applied to estimate the effective properties of random 2D composites. 
First, deterministic boundary value problems are solved for all loca¬ 
tions of inclusions, i.e., for all events of the considered probabilistic 
space C by the generalized method of Schwarz. Second, the effec¬ 
tive properties are calculated in analytical form and averaged over C. 
This method is related to the classic method based on the average 
probabilistic values involving the ?i-point correlation functions. How¬ 
ever, we avoid computation of the correlation functions and compute 
their weighted moments of high orders by an indirect method which 
does not address to the correlation functions. The effective proper¬ 
ties are exactly expressed through these moments. It is proved that 
the generalized method of Schwarz converges for an arbitrary multi¬ 
ply connected doubly periodic domain and for an arbitrary contrast 
parameter. The proposed method yields effective in symbolic-numeric 
computations formulae of high order in concentration. In particular, 
the Torquato-Milton parameter £1 is exactly written for circular in¬ 
clusions. 


1 Introduction 

The effective properties of random heterogeneous materials and methods of 
their computation are of considerable interest n. a, a, ns. P3, m, 
Hu. m- Randomness in such problems reveals through random tensor- 
functions locally describing the physical properties of medium. Despite the 
considerable progress made in the theory of disordered media, the main tool 
for studying such systems remains numerical simulations. Frequently, it is 
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just asserted that it is impossible to get general analytical formulae for the 
effective properties except dilute composites when interactions among inclu¬ 
sions are neglected and except regular composites having simple geometric 
structures. This opinion sustained by unlimited belief in numerics has to 
be questioned by the recent pure mathematical investigations devoted to 
explicit solution to the Riemann-Hilbert problem for multiply connected do¬ 
mains [24], [32], [35! and by significant progress in symbolic computations 
[9], [10]J. In the present paper, we demonstrate that the theoretical results 
m, E2). [35] can be effectively implemented in symbolic form that yields 
analytical formulae for random composites. 

For conductivity problems governed by Laplace’s equation, the local con¬ 
ductivity tensor A(x) can be considered as a random function of the spa¬ 
tial variable x = (xi,X2,x 3 ). In the present paper, we restrict ourselves to 
two-phase composites with non-overlapping inclusions when a collection of 
particles with fixed shapes and sizes is embedded in matrix. More precisely, 
let be given a set of hard particles {V k , k = 1,2,...} where each T\. has 
a fixed geometry. Let all the particles are randomly located in the space 
and each particle T> k occupies a domain D k without deformations. Thus, 
the deterministic elements D k are introduced independently but the joint set 
{Di, LL, • • •} randomly. The diversity of random locations is expressed by 
joint probabilistic distributions of the non-overlapping domains D k . 

Various methods were applied to such random composites. The most 
known theoretical approach is based on the n-point correlation functions 
presented by Torquato 021 and summarized below in Secj3] This theory 
yields analytic formulae and bounds for the effective properties. The main 
shortage of this theory is computational difficulties to calculate the n-point 
correlation functions for large n. 

Let a distribution of inclusions be fixed. Then, the effective conductivity 
tensor A e is uniquely determined through T> k and their physical properties. In 
order to estimate A e , various statistical approaches were applied to construct 
the representative volume element (RVE) [T6]. These approaches are based 
on the straight forward computations of the effective properties for various 
locations of inclusions and on the statistical investigations of the obtained 
numerical results. The main lack of this physical theory is that the notation 
of the RVE is not correctly defined and the results have not precise meaning. 
This theory is rather a statistical analysis of the computational experiments 
and physical measures. 

The main purpose of this paper is to work out constructive analytical- 
numerical methods to calculate the effective constants of composites with 
non-overlapping inclusions and to develop the corresponding RVE theory. 
We apply the straight forward approach to estimate the effective conductivity 
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of random 2D composites. First, deterministic boundary value problems are 
solved for all locations of inclusions, i.e., for all events in the considered 
probabilistic space C by the generalized method of Schwarz (GMS). This 
method was proposed by Mikhlin 125 as a generalization of the classical 
alternating method to multiply connected domains; convergence of the GMS 
was established in [25], [35]. The GMS is referred to decomposition methods 
[40 ] frequently used in numerical computations and realized in the form of 
alternating methods. A method of integral equations closely related to the 
GMS was developed in |8j. 

After solution to the deterministic problem by the GMS the ensemble av¬ 
erage for the effective conductivity (mathematical expectation) is calculated 
in accordance with the given distribution. This method is related to the 
classic method [42J based on the average probabilistic values involving the 
n-point correlation functions. In our approach following the GMS, we avoid 
computation of the correlation functions and compute their weighted mo¬ 
ments. The effective properties are exactly written through these moments. 
Analytical formulae and numerical simulations demonstrate advantages of 
our approach in the case of non-overlapping inclusions. For instance for 
circular inclusions, our method yields effective in computations formulae of 
order 0(is 2 °) in concentration v (exactly written in Sec 14. 5 1 up to 0(v 8 )). 
Constructive formulae of the classic approach 02! are deduced only up to 
0(u 5 ) for not high contrast parameters. In particular, the Torquato-Milton 
parameter (j [22], [52] is exactly written for circular inclusions. 

Stochastic 2D problems are posed and solved in doubly periodic statement 
in the plane. Theoretically, doubly periodic problems constitute the special 
class of problems in the plane with infinite number of inclusions 125 - How¬ 
ever, the number of inclusion per a periodicity cell, N, is arbitrary taken, and 
the final formulae for the effective tensor contains N in symbolic form. Sim¬ 
ilar non-periodic statements can be used following [34]. Such an approach 
explains, for instance, the divergent sum (integral) arising in applications of 
self consistent methods when the divergent sum does really diverge for some 
distributions. This implies that the corresponding composites cannot be ho¬ 
mogenised over the whole plane even if the concentration is properly defined. 
The homogenization theory (see [41] and works cited therein) justifies exis¬ 
tence of the effective properties for statistically homogeneous random fields 
which constitute a subclass of heterogeneous fields discussed in 021. 122], 
[34]. The divergence and other similar effects do not arise for doubly peri¬ 
odic composites when homogenization is always valid. A periodicity cell can 
be treated as a representative cell and vice versa. From practical point of 
view, each sample is finite, hence, it can be considered as a representative 
cell with many inclusions. Application of the new RYE theory (see the next 
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paragraph) can reduce the number of inclusions per periodicity cell but not 
always significantly. Because reduction to a small number of inclusions per 
cell can distort the effective properties of the original sample. It was noted 
in [6J that periodic arrange of inclusions decreases the effective conductivity 
when conductivity of inclusions is greater than conductivity of host. More¬ 
over, the theory of elliptic functions can be applied to constructive symbolic 
computations for periodic problems that yields exact and approximate ana¬ 
lytical formulae for the effective tensor. 

The proposed method leads to construction of the rigorous RVE theory 
for non-overlapping inclusions [30]. The effective tensor can be written in 
the form of expansion in the moments of the correlation functions which can 
be considered as ’’basic elements” depending only on locations of inclusions. 
The RVE is defined as the minimal size periodicity cell corresponding to the 
set of basic elements calculated for the composite. A simple fast algorithm to 
determine the representative cell for a given composite is based on reduction 
of the number of inclusions per periodicity cell having the same basic elements 
as the given composite. It follows from simulations i. H that for the 
uniform non-overlapping distribution of disks in order to reach high accuracy 
in the effective conductivity one needs to solve the corresponding boundary 
value problems at least for 64 inclusions per cell repeated at least 1500 times. 
These technical parameters essentially exceed possibilities of the traditional 
statistical approach to the RVE. 

Here, we are restricted by composites with non-overlapping inclusions. 
Other types of composites and porous media, for instance when A(x) obeys 
the Gaussian distribution, the correlation functions and reconstructions of 
media give excellent practical results (see |1] and works cited therein). 

The goal of this paper is to show that the results obtained by the GMS 
a. ra. m-m, m can be considered as a constructive approach to re¬ 
solve some problems of the classic theory of random composites ra for non¬ 
overlapping inclusions. The GMS allows to overcome myths frequently met 
in the theory of random composites: 

1. It is asserted that only application of the correlation functions is the 
proper method to treat random composites. Hence, restrictions in this field 
are related to computational difficulties arising in numerical computations 
of the correlation functions. As it is noted above, we obviate the need for 
the correlation functions calculating the mathematical expectation of the 
effective tensor after complete solution to the deterministic problems. 

2. It is thought that the contrast parameter expansions are valid only 
for sufficiently small contrast parameters. However, it was first justified in 
[23] that such expansions for plane conductivity problems are valid for all 
possible values of the contrast parameters. 
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3. Self consistent methods for random composites are considered as a 
simple and constructive alternative to the method of correlation functions. 
However, applications of the GMS demonstrates essential limitations of these 
methods [34]. More precisely, the self consistent methods do not capture 
interactions among inclusions. Even solution to the n-particle problem in 
the infinite plane can be applied only to diluted distributions of the clusters 
consisting of n inclusions. 

This paper is the Erst part of the work devoted to the GMS described 
in general form in Sec l4.ll The special attention is paid to conductivity 
problems. Integral equations corresponding to the GMS are described in 
Sec l4.21 These equations can be solved by different iterative schemes. The 
first scheme is based on the contrast expansions (the classic approach is in 
Sec l3.2l and the GMS form in Sec 14.31) . The second method uses the cluster 
expansions (the classic approach in Sec 13.II and the GMS form is in Sec 14. 41) . 
Applications of the GMS to elasticity problems will be described in a separate 
paper. 

2 Deterministic and stochastic approaches 

Let uj i and 0 J 2 be the fundamental pair of periods on the complex plane C 
such that oj\ > 0 and Im U 2 > 0 where Im stands for the imaginary part. The 
fundamental parallelogram Q is defined by the vertices ±Wi/2 and ±cu 2 /2. 
Without loss of generality the area of Q can be normalized to one, hence, 


( 2 . 1 ) 


uqlm uj 2 = 1. 


The points m\ 0 J\ + rri 2 UJ 2 {mi, m 2 G Z) generate a doubly periodic lattice Q 
(torus) where Z stands for the set of integer numbers. 

Consider N non-overlapping simply connected domains D k in the cell 
Q with Lyapuniv’s boundaries L k and the multiply connected domain D = 
Q\ (D k U L k ), the complement of all the closures of D k to Q. Each 
curve L k leaves D k to the left. Let a point a k is arbitrary fixed in D k (,k = 


1,2 


We study conductivity of the doubly periodic composite when the host 
D + miUi + m 2 0 J 2 and the inclusions D k + m l oj\ + m 2 uj 2 are occupied by con¬ 
ducting materials. Without loss of generality the conductivity of matrix can 
be normalized to one. The respective conductivity of inclusions is denoted 
by A. Introduce the local conductivity as the function 



z e D, 

z E D k , k = 1, 2,..., N. 


( 2 . 2 ) 
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where z — x\ + ix 2 denote a complex variable. The potentials u(z) and u k (z) 
are harmonic in D and D k (k — 1,2, ..., N ) and continuously differentiable in 
closures of the considered domains. The conjugation conditions are fulfilled 
on the contact interface 

u = u k , on L k , k — 1,2, ..., N, (2.3) 

on on 

where d/dn denotes the outward normal derivative to L k . The external held 
is modelled by the quasi-periodicity conditions 

u(z + ui) = u(z ) + £i, u(z + u> 2 ) = u(z) + £ 2 , (2.4) 

where £^ (£ = 1, 2) are constants. The external held applied to the considered 
doubly periodic composites is given by the vector E 0 = — (£i,£ 2 )- 

The problem (I2.3ll -' |2.4ll can be stated in the probabilistic context with 
randomly located inclusions. For definiteness, consider a set of the domains 
D k as a realization of the random locations of a collection of particles V k with 
hxed shapes and sizes. It corresponds to location of hard particles in a vessel. 
Following the definition (3.4) from [42] we introduce the doubly periodic 
specihc probability density P(a) associated with finding a configuration of 
inclusions with position a := (di,a 2 ,... ,ajy). The point a k determines the 
position of the particle T> k in the cell Q. The set of all the configurations is 
denoted by C. Let da*, denote a small area element about the point a k and 
da = daida 2 ... da at. Then P(a)da is equal to the probability of finding a± 
in dai, a 2 in da 2 , ..., ajv in dajv. 

The ensemble average is introduced as the expectation in a probabilistic 
space 

/ = f /(a)P(a)da. (2.5) 

The effective conductivity tensor of deterministic composites 

Ae = ( ) (2.6) 

\ / '21 A 22 J 

can be found in terms of the solutions of two independent problems (12. 3 p - 
(12.4p . For instance, one can take £1 = uji, £ 2 = Re co 2 and £1 = 0, £ 2 = Im uj 2 . 
In the first case [7], 


N 

An = 1 + 2 p 

k =1 


duk 

dx\ 


dxidx 2 . 


(2.7) 
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Here, p denotes the contrast parameter introduced by Bergman 


P 


A- 1 

ITT 


( 2 . 8 ) 


The effective conductivity tensor A e of stochastic composites is deter¬ 
mined by application of the operator (12.5)1 to the deterministic values (12.6p . 
Let the probabilistic distribution of the domains D *. (k = 1,2,..., N) in 
Q be such that the corresponding two-dimensional composite is isotropic 
in macroscale, i.e., the expected effective conductivity of the composite is 
expressed by a scalar. It can be calculated as follows 


A e — An — / An(a)P(a)da, (2-9) 

where the deterministic value An (a) is calculated by (12.711 . 


3 Classic approach 

As it is stressed in Introduction we apply the straight forward approach to 
determine A e . First, the deterministic problem (I2.3j) (I2.4j) is solved for every 
configuration a £ C by the GMS. Further, the probabilistic average (12.911 is 
calculated. In order to compare the results obtained by this method and by 
others we follow the classic approach to random composites described and 
summarized by Torquato [32]. We follow the book [42] and do not present 
formulae from earlier classic works (see references and historical notes therein 
and also in 1221) in order to avoid multiple notations. We stress that only 
two-dimensional two-phase dispersed composites with one phase embedded 
in the connected matrix are considered in the framework of the general theory. 

In the book [42], the problem (12.311 -- (12. 1 D is written in terms of the electric 
conduction by use of the local electric current J and the intensity field E. 
The potential u satisfies equation 

E = — Vm. (3.1) 

The linear constitutive relation is fulfilled in D and Dk 

J (z) = A(z)E(z). (3.2) 

Torquato 021 and others distinguish two main methods: cluster expansion 
and contrast expansion discussed below. 
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3.1 Cluster expansions 

This method is presented as a correction of the single-inclusion approach 
by taking into account interactions between pairs of particle, triplets and so 
forth. It is worth noting that only infinite number of inclusions on plane gives 
a correct result [42], [33]- Study of the finite number of inclusions on plane 
yields the effective properties only of dilute composites as it is demonstrated 
in p>41| . 

The intensity field E is considered as a function of z with the configuration 
parameter a and denoted by E(z;a). Let the constant external field Eo is 
fixed. Torquato describes the following cluster expansion (see (19.8) from 

021 ) 


N N N 

E(z; a) = Eo + ^ Mi(z; k) ■ Eo + ^ M 2 (z; k , k\ ) • Eo + • • • . (3.3) 

k= 1 k =1 k\y^k 

Here, Mj(z; fc) is the single-inclusion operator which accounts for the first 
order interactions over and above E 0 and inclusions periodically equivalent to 
the fc-th inclusion, i.e., to D k + miUJi + 77220)2 for all integer m\ and m 2 • The 
operator Mi(z;fc) can be derived by the following two methods. The first 
method corresponds to the Maxwell approach applied to dilute composites. 
It consists of two steps. At the first step, a boundary value problem for one 
inclusion D k in the infinite plane is solved and the single-inclusion operator 
K i(z]k) in the plane is constructed (see for instance, formulae (17.4) and 
(19.6) from [32]). Next, Mi(z; k ) is constructed via Ki(z; k ) by a periodicity 
operator. The periodicity transformation V : Ki(z;fc) t—)■ Mi(z; fc) can be 
easily performed since K^z; aj) is expressed via the dipole tensor. The 
periodicity operator transforms rational functions to elliptic functions [2); 
for instance, V : z~ x t-)- £(z) and V : z~ 2 t-)- p(z). The second method to 
construct Mi(z; k ) is based on the same scheme realized in the torus topology 

PEj- 

The operator K^z; k ), hence M^z; k ), is written in closed form for disk, 
ellipse and many other domains (see [32] an d works cited therein). It can be 
constructed for any shape D k in terms of the conformal mapping of D k onto 
the unit disk. In general, K p (z; k, ki,... , k p ) and M p (z; k, k\,... , k p ) are the 
p-inclusion operator which accounts the pth order interactions. 

The effective conductivity can be found through the polarization field 
defined by (19.3)-(19.4) from [ 32 ] 


P(s) = [A( 2 ) - 1]E(2) 


(A — l)E(z), z G D , 

0, z G D k , k = 1,2,..., N. 


(3.4) 


Following the classic approach [32] we introduce the ’’double” average oper¬ 
ator over the unit cell Q and over the configuration space C 


{g)= / g(x i,x 2 ;a)dadxida;2. 


Q Jc 


The effective tensor (12.61) can be defined by equation 


(P) = (A e -I)-(E), 


(3.5) 


(3.6) 


where I is the unit tensor. Having used (13.61) Torquato (see (19.33) from [32] 
and earlier investigations cited therein) deduced formula for macroscopically 
isotropic composites 


A e 



114(a) da, 


(3.7) 


where W n { a) is a complicated functional of the n- inclusion cluster operators 
M n (z; k,ki,... i k n ) and the probability density P(a). It is worth noting that 
W n (a) is derived through K. n (z‘, k,k\,..., k n ) in [32], Actually, our periodic 
approach yields the same result by rearrangement of the terms in each W n { a). 

Let \Dk\ denote the area of Dk■ The concentration of inclusions is defined 
in the following way 

N 

u = J2\ D kl (3-8) 

k =i 

It is noted in [32] that (13.71) yields the volume-fraction expansion (see (19.34) 
from [32] ) 

oo 

\ e =l + J2 B n" n , (3-9) 

n =1 

where the coefficients B n are multidimensional integrals on M p (z; k, Aq,..., k p ) 
for p = 1,2,... ,n. Direct applications of the correlation functions yield an¬ 
alytical formulae and numerical results for B n n = 2,3,4 (see 02 and a 
discussion in [9]). 


3.2 Contrast expansions 

Exact contrast expansions can be presented by formula (20.1) from [32J which 
after some modifications becomes 

OO 

\ e =l + J2 C nP n , (3-10) 

71=1 
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where p is the contrast parameter 02.81) . Many authors have been thinking 
that the series (13. 10p converges only for sufficiently small p. Though conver¬ 
gence of the series for local fields had been proved in 12a. 12a for circular 
inclusions and in [35] (see Sec.4.9.2) for general shapes for all admissible 
\p\ < 1. This of course implies the convergence of 03.101) . 

In order to deduce the general representation 03. lOj) and estimate the first 
coefficients C n (n = 1,2, 3,4) Torquato [42] derives an integral equation on 
E(z) which can be written in the form (compare to equation (20.17) from 
021 for the cavity intensity field) 

E (z) = E 0 + f d z' G (z, z') • [A (z) - 1]E(0, (3.11) 

Jq 

where G(z, z 1 ) denotes the periodic Green’s function. Substitution of 02.2j) 
into 03.lip yields 


E(z) = E 0 + p(X + 1) Y" / dz' G(z, z') ■ E(/). (3.12) 

fc=i Jd * 

The method of successive approximations applied to 03.121) yields a series on 
the contrast parameter p which can be shortly written in the form 

OO 

E(z) = p n H n ■ E 0 , (3.13) 

n= 0 

where the operator H is dehned by the right hand part of 03. 121) . Use of 
03.4p yields an analogous series for P(z). The tensor A e satisfies equation 
03.61) which can be written in the form (20.37) from [42]• The latter for¬ 
mula for macroscopically isotropic media becomes (20.57) from [42]. In our 
designations it reads as follows 

OO 

pV(A e - l) -1 (A e + 1 ) = vp-^2 A^p n . (3.14) 

n =3 

The n-point coefficients An ^ are exactly expressed in terms of the p-point 
correlation functions (p — 2,3, , n ) (see formulae (20.38)-(20.41) and (20.59)- 
(20.60) from [12jJ)- Formula (j3. 14j) can be written in the form (13.101) . 

The coefficients An > were calculated in [12] as sums of multiple integrals 
containing the correlation functions which are presented also through multi¬ 
ple integrals. Calculation of these integral is the main difficulty for numerical 
application of (j3 . 14j) . 
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4 The generalized method of Schwarz 

4.1 General 


The generalized method of Schwarz (GMS) was proposed by Mikhlin [21] as a 
generalization of the classical alternating method of Schwarz for finitely con¬ 
nected domains. The GMS is based on the decomposition of the considered 
domain with complex geometry onto simple domains and subsequent solu¬ 
tion to boundary value problems for simple domains. In our case, we have TV 
simple domains, the cell Q with one inclusion D k (k = 1, 2,..., TV). In each 
step of the algorithm the boundary conditions for the fcth simple domain are 
corrected by influence of the mth simple domains (m ^ k) computed at the 
precedent steps. 

We now shortly present the GMS, first for the conductivity problem (12.31) - 
(12. 4p and further for general linear problems. The ideal contact condition 
(12.31) can be rewritten in the form 


u o = u k -u 


ext 


duo 

dn 


d(u k - u 


ext\ 


dn 


+p(A+l) 


du k 

dn 


on Li 


k = 1, 2,..., TV, 


(4.1) 

where u = w 0 + u ext is the decomposition of the potential u onto the regular 
periodic part u 0 and the external field u ext corresponding to the quasi peri¬ 
odicity conditions (12.4[) . In the case of the space problem the external field 
u ext has a singularity at infinity. 

Introduce the function 


f(t) = (A + l)^(t), t e L k (k = 1, 2,..., TV) 


and the domains D + := U k=1 D k , D~ := D. It is convenient to represent the 
harmonic functions in different domains as one piece-wise harmonic function 
uq(z) introduced above in D = D~ and equal to u k (z) — u ext (z) in D k U L k 
(k — 1,2,, TV). Then, (14.11) can be considered as the jump problem 


(itI ' (itI 

<=u- Ql -^L = -£ + p f( t ) on L = U k=1 L k , (4.2) 

where Uq and Mq correspond to the limit values of u(z) when z tends to t 
from D + and D~, respectively. The jump problem (14.21) has a unique solution 
expressed in terms of the simple layer potential P for the curve L in the cell 
Q in the torus topology. We have 


u 0 (z)=p(Pf)(z ), ZED+UD-. (4.3) 
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The operator P is considered as an operator in an appropriate f un ctional 
space HD, HD, Id, M- In the classic statement, Q should be replaced 
by the whole space. The operator P is decomposed onto the simple layer 
potentials Pk along curves L k (k — 1, 2,..., N) as P = J2k=i -Hfc- Here, 
the multiplier yd_- j s introduced for short form of the equations below. Then 
(14.31) implies that 

N f A,,+ \ 

U k (z) = pJ2 { p m-Q^ J (z) + U ext (z), zeD k (k = l,2,...,N), (4.4) 

m= 1 ' ' 


and 

JT / rh,+ \ 

u(z) f - p m-Q^ J (z) + n ext (z), z E D. (4.5) 

m= 1 ' ' 

Equations (14. 4 1) can be considered as a system of integral equations on the 
potentials u k {z) in the inclusions D k (k = 1,2If it is solved, the 
potential u(z) is calculated in D by (14.5[) . 

Remark 4.1. Equations (14.4j) correspond to the GMS in Mikhlin’s form ED 
when the method of successive approximations can diverge for closely spaced 
inclusions. The following slight modification yields a convergent algorithm 
for any location of inclusions and every contrast parameter (|p| < 1) 



U) - ( P m ) 


(w) +u ext (z)+ 


z E D k L) L k (k — 1, 2,..., N), 


(4.6) 


where w is a fixed point in D (see explanations in Sec.4 of [32]). 

Equations (I4.3j) — (14.61) are written for potentials. Their differentiation 
yields equations on the intensity field 


N 

Efc(z) = p ^^(Q m E m )(z) + E 0 (z), zED k UL k {k — 1, 2 ,..., N), (4.7) 

m =1 


where E k (z) = E(^) in D k and Q m are appropriate operators (similar to 
double layer potentials). Equation (14.71) is similar to (13. lip . 

This GMS for harmonic functions can be extended to general problems 
governed by linear equations describing electrical and heat conduction, flow 
in porous media, viscous flow, elasticity and coupled fields in K d (d = 2, 3). 

Let tensor fields U(z) in D and Ufc(z) in D k satisfy the contact conditions 


T • U = Tfc • Ufc, on L k (k = 1, 2 ,..., N), (4.8) 
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where T and T& are boundary operators involving physical parameters of 
materials occupying the domains D and D k , respectively. Let T fc is presented 
in the form Tfc = T— p k , where p k denotes a contrast parameter tensor. Then 
(14.8)1 becomes 

U = U fe — T -1 • p k • Ufc, on L k {k = 1,2,... ,7V), (4.9) 

Let Pfc denote the simple layer potential corresponding to the considered 
linear equation for the domains D k and D k separated by L k . Application of 
Yhk=i Pfc to (S3D yields 

N 

Ufc ( 2 ) = J2 ( p m • T- 1 • p m ■ U m )(z) + U ext (z), z E D k UL k (k = 1,2,..., N). 

m= 1 

(4.10) 

There are two different methods to solve equations (14 .1 01) . The hrst 
method is based on the direct iterations and corresponds to the contrast 
expansion in p (compare to Sec l3.2l) . The second method is based on implicit 
iterations applied to the same equations (14.101) written in the form 

U t (z) - (P t . T- 1 ■ p k . U*)(z) = ELt( P ". ' T_1 ■ Pra ■ UJ W + U 

z £ Dk U L k {k = 1 , 2 ,..., N). 

(4.H) 

One-inclusion problem for D k (k = 1,2,, N ) is solved at each step of 
iterations. This scheme corresponds to the cluster expansions discussed in 
Se d3.ll 

Remark 4.2. Equations (14.31) and (14.111) are not reduced to the classic integral 
equations constructed in the framework of the potentials theory [T9] . It is 
another type of equations. 

Remark 4.3. Equations (I4.3f) — (14.51) and (14.101) - (14.11)) can be derived in the 
limit cases of inclusions (soft and hard inclusions in terminology [38]) by 
introduction of hctive potentials [35]. 

4.2 Integral equations 

In the present section, we discuss the GMS for the double periodic problem 
(I2.3l) - (j2.4j) . Following [36] and [7] we introduce the complex potentials cp(z) 
and 0fc(z) in such a way that 

u (z) = Re <f>(z), u k (z) = —y—-Re (z) . (4.12) 

A + 1 
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Two real equations (12. 3 p can be written as one R- linear condition (see for 
details [2Tj, [36] and [7]) 


(j) ( t ) = (p k (t) - p (pk (t), te L k , k = 1, 2,..., N, (4.13) 

where the bar denotes the complex conjugation. The unknown functions 
(p (z) and (pk (z) are analytic in D and D k , respectively, and continuously 
differentiable in the closures of the domains considered. It follows from (I2.4[) 
that the function <p(z) is quasi periodic |28j 

<t>{z + Uj) - <j>(z) = + idj (j = 1,2), (4.14) 

where i stands for the imaginary unit, the constants £i and £2 are taken 
from (I2.4[) . d\ and d 2 are undetermined real constants which should be found 
during solution to the problem. 

The following formula was proved in [35] 

= ~n 2 {t) (p' k (t), t G L k . (4.15) 

Here, the differential operator is understood as ^ = t' s j-, where t = t(s ) is 
a parametric equation of L k and 4>' k (t) is the boundary value of the complex 
derivative — i-£^- Differentiation of (j4.13j) and application of (I4.15P 

yields 

-0 (t) = (t) + pn 2 k (t) ip k (t), t e L k , k = 1,2,..., N, (4.16) 

where ip(z) = 4>'(z), 'ipk(z) = (p' k ( z ) and n k (t) stands for the unit normal 
outward vector to L k expressed in terms of the complex values. The function 
ip(z) is doubly periodic. The M- linear problem (14.161) holds also in the limit 
cases for perfectly conducting inclusions (p = 1) and for isolators (p = —1). 

The M- linear problem (14.131) - (II. 1 II) can be reduced to a system of integral 
equation in the following way. The cell Q in the torus topology is divided 
onto two domains D + = U k=1 D k (not connected) and D~ = D (multiply 
connected). Let L = U k=1 L k denote the boundary of D + and p(t) be a 
Holder continuous function on L. Let £( 2 ) denote the Weierstrass ((-function 
for which |2J 

C(z + Uj) -CO) = Sj (j = l,2), (4.17) 

where Sj = 2( (^-). Following [43], [37] introduce the Eisenstein function 

Ei(z)= y ---, (4.18) 

^ v z + miOJi + m 2 uJ2 

mi 
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where the Eisenstein summation is used [33]. The Weierstrass and Eisenstein 
functions are related by formula Ei(z) = ((z) — S 2 z, where S 2 = Jj-C (it) • It 
follows from Legendre’s identity 8\uj 2 — S 2 oji = 2m (see [2], PJ) and (14471) 
that the jumps of E\{z) have the form 

27 

Ei(z + uq) — Ei(z) — 0, E 1 (z + u 2 ) - Ei(z) = -. (4.19) 

OJi 

The Cauchy-type integral on torus was introduced with the ((-function in 
the kernel [3]. We introduce it in the slight different form using the Eisenstein 
function 

<&(z) =- [ n(t)E\(t — z)dt + Cz + Co. z G D ± , (4.20) 

2 vrz J L 

where fi(t ) is a Holder continuous function on L, C and Co are complex 
constants. The function $( 2 ) is analytic in and quasi-periodic, i.e., it 
has constants jumps per periodicity cell 

&(z + uji) — $(z) = Coji, $>{z + cu 2 ) — $( 2 ) = — f p(t)dt + Cu 2 . (4.21) 

wi Jl 

The following Sochocki (Sokhotski-Plemelj) formulae take place on torus [3] 

® ± ( t ) = ± \v( t ) + ^jJri T )Ei{T-t)dT, teL, (4.22) 

where < h ± (t) denote the limit values of < h(^) on L when z G D ± tends to 
t G L, respectively. Formulae (14.221) imply that the function 4*(^) satisfies 
the jump condition 

4> + (t) - $"(t) = t G L. (4.23) 

Equation (I4.23P can be considered as a C-problcm problem in a class of 
quasi-periodic functions with the given jump p{t). It follows from [3J that 
the general solution of (j4.23j) up to an additive arbitrary constant has the 
form (14.201) . 

Consider (14.13j) as the problem (14.23[) with fj,(t) = p 4> k (t) on L, <I>(z) = 
4>k{z) in D k C D + and $(z) = c f>{z ) in D = D~. Then, (14.20)1 can be written 
in D + as follows 

N 1 r _ 

4>k{z) = / ^> m (t)E 1 (t-z)dt + Cz + C 0 , z G D k (. k = 1,2,..., N). 

ro —1 J Lm 
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(I4.20p in D becomes 


<t>(z) = ^7 f (j) m (t)Ei(t - z)dt + Cz + C 0 , z e D. (4.25) 

m =1 

Equations fj4.24h - fl4.25h are deduced from the R-linear problem (I2.3p ~ fl2.4p . 
Using the same arguments one can derive analogous equations from the prob¬ 
lem (14. 161) on the complex flux 





n 2 (t ) 'ijj m (t)Ei(t—z)dt+C, z e D k (k 


1 , 2 , 


^(z) 



n 2 (t) 


z)dt + C, z G D. 


(4.26) 

(4.27) 


It is also possible to obtain these equations by differentiation of (j4.24p ~ fj4.25p 
using integration by parts and formula (I4.15p . 

The integral equations (14.261) are considered as an equation in the follow¬ 
ing Banach space. First, the space (U° ,Q 0 Q f functions Holder continuous on 
L is considered (0 < a < 1). This is a Banach space endowed with the norm 


\UJ 


sup I UJ (t) | + sup 

t£L t\ : 2&L 0^*2 


u(tl) - <4*2)1 
\tl ~t 2 \ a 


(4.28) 


Consider the closed subspace A C C^ 0,a ' > consisting of functions analytically 
continued into all D k . Therefore, every element of A is a function analytic in 
the domain U k=1 D k and Holder continuous in its closure. It is worth noting 
that the domain U k=1 D k is not connected and convergence in A implies the 
uniform convergence of functions in U k=1 (D k U L k ). 

Equations 04.261) are written in the domains D k and can be continued to 
the boundary by use of (14.221) 




-p 2^ra= 1 


— jip 1 
2 IL m\ 


(t) i>. m(t ) + Jr n 2 m (t) -z)dt + C, 


t E L k (k — 1,2,..., N). 

(4.29) 

Therefore, the integral equations in the space A have the form 04.26p . (14.291) . 
It is proved in Appendix that equations (I4.26p . (I4.29P have a unique solution. 
Solution to these equations corresponds to the GMS which can be realized 
in different forms as it is noted in Sec 14.II 

The problem (14.13p (14. 14p for the complex potentials has a unique solu¬ 
tion up to an arbitrary additive constant Cq. The constants £i + id\, £2 + hi 2 
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and Co disappear in the problem (I4.16p after differentiation. The structure 
of the general solution of (14.16jl has the form 'ip(z) = £1 + £ 2 pA 2 ' 1 (z) 

(see [28]), where the real constants £1 and £2 correspond to the external field 
(jumps of u(z) per a cell). Hence, the M-linear problem (14.16]) has two M- 
linear independent solutions. From the other side, these two independent 
solutions can be constructed by (j4.26|) (j4.27[) taking, for instance C = 1 and 
C = i. Therefore, the complex constant C can be expressed through the real 
constants £1 and £ 2 - In order to calculate the effective properties of rnacro- 
scopically isotropic composites it is sufficient to fix £ 1 , £2 and fold C. For 
simplicity, we put 

£1 = kh, £2 = Re 0 J 2 - (4.30) 

Then the external field corresponds to the complex potential £>( ext )[ z ) = z 
and to the complex flux ■£?( ext )(~) — 1 . j n the vector form, this external flux 
becomes (—1,0). It follows from the first relation (14.21 p that 


C 


i+hi, 

Ul 


(4.31) 


where d\ 
calculate 

d\ — Im 


= Im [<j)(z + Wi) — (j)(z )]. Put z = — \{oj\ + UJ 2 ) in this relation and 


r 5(0)1—w 2 ) 
' — +^ 2 ) 


if;(t)dt = — 


i-Reoj 2 ) Q u 
i+RiJ2) &X‘2 


ilm ujo ', 

xi --— ] dxi. (4.32) 


The integral from the left hand part of (14.32(1 expresses the total flux through 
the low side of the parallelogram Q parallel to the :r r - axis. It must be equal 
to zero because of the periodicity and that the external flux (—1, 0) vanishes 
in the .xy-direction. Therefore, d\ = 0, hence C — 1 by (14.311) . 

The effective conductivity tensor of deterministic composites can be cal¬ 
culated through the complex potentials in the inclusions (see [27], [36] and 

m) 

N 

An—iAi2 = l + 2 ipk{xi + ix2) dxidx 2 , (4.33) 

k =1 ■ Dk 

where ipk(z) (k — 1, 2,..., N) satisfy (14.26(1 — (]4.27() . (I4.29|) with C — 1. The 
value (I4.33|) depends on the locations of inclusions, i.e., on a 6 C. The effec¬ 
tive conductivity of macroscopically isotropic random composites is defined 
through the operator (]2,5jl applied to (14.331) 


A e 



(4.34) 


where the macroscopic isotropy implies that A 12 vanishes and An = A 22 - 
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4.3 Contrast expansions for the GMS 


It is proved in Appendix that the unique solution of (14.261) . (14.29)) can be 
found by a method of successive approximations uniformly convergent in 
U^ =1 (Dk U Lj.). Let T(z) = for z G D k ALk and A denote the operator 

from the right hand part of ( 14.24ft . (14.291) . Then, equations (14.241) . ()4.29j) can 
be shortly written as the following equation in the space A 


T = pAT + 1. 

Its unique solution is given by the series 


T = ^p n A n 1. 


(4.35) 


(4.36) 


n=0 


The method of successive approximations can be also presented in the 
form of the following iterative scheme 

'!>k\ z ) = h _ 

^t\ z ) = - z)At + 1 , p = 1 , 2 ,..., 


zeD k (A: = 1,2,..., TV). 

(4.37) 

When the functions if)k(z) (k — 1, 2,..., N) are found, the effective conduc¬ 
tivity is calculated by (14.330 . 

In order to compare equations ()3.12j) and (14.241) . (14.291) we introduce the 
complex intensity held E(z) = E\ (z) — iE 2 (z ) isomorphic to the vector held 
E(z) = (Ei(z), E 2 (z)). Then, the vector equation (13.121) can be written 
as a complex equation on E(z). For dehniteness, we consider the complex 
potential V'fc(^) in D k . It follows from (13.11) . (14.121) and "0fc(^) = (p'k( z ) that 


E(z) = 


dui 


.dm 




(4.38) 


dxi dx 2 A +1 
Therefore, the complex equation on E(z) in D k can be written as the equa¬ 
tion (14.24)1 on ijjk(z). We do not directly prove that equations (13.121) and 
(14.35)1 are the same. But it is easy to show that the forms of their solutions 
(13.13)1 and ()4.36jl coincide. Since the solutions E and T coincide (more pre¬ 
cisely, isomorphic) for E 0 = (—1,0), the power series (13.13ft and (14.361) in p 
can be considered as the presentations of the same function for sufficiently 
small \p\. Therefore, the coefficients of (I3.13P and (14.361) coincide. In partic¬ 
ular, this fact implies that substitution of ()4.36jl into ()4.33jl yields a series 
which coincides with the series (13.101) obtained by the classic approach. It 
follows from ra, m, |.35j that the contrast expansion (14.36)1 converges for 
all admissible |p| < 1. 
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4.4 Cluster expansions for the GMS 

The GMS can be presented by other iterative scheme corresponding to the 
classic cluster expansions summarized in Sec 13. 11 First, we rewrite (I4.24jl in 
the form 

M z ) + 3^ f Lk n l(*) ~z)dt = 


-pT,m^k^riI Lni n m(t) — z)dt + 1, z e D k (k — 1,2,... ,n). 

(4.39) 

Following Sec l4.2l for every fixed m = 1,2,..., iV introduce the Banach 
space A m of functions analytic in D m and Hblder continuous in D m U L m . 
Let f m (z ) G A m - For each fixed k — 1,2,..., iV introduce the operators 


( Pkmfm){z) = -— rf n {t) f m (t)E\(t - z)dt, z G D k (m = 1 , 2,..., n). 

Am l T 

u J-'m 

(4.40) 

One can check that the operator P kk is singular and P km {jn ^ k) are compact 
in Ak- Equations (|4.39[i can be shortly written in the form 


(/ + pPkk)^k = ~(> y Pkmi’m + 1, k = 1,2,... ,n. (4.41) 

m^k 

The zero-th approximation in the concentration v for (14.39(1 can be written 
as the following integral equation 


t \ z ) + ■£- f nl(t) ip^\t)Ei(t - z)dt = 1, z G D k . (4.42) 

Equation (14.421) corresponds to the M- linear problem (I4.16|) for one inclusion 
D k in the cell Q. Solution to the one-inclusion problem (|4.42[) defines the 
inverse operator (/ + pPkk)~ l bounded in A k - Then, (14.411) can be written in 
the equivalent form 

ipk = -p ^ (I + pPkk)~ 1 Pkm^m + (/ + pPkk)~ l 1, k = 1, 2,. .., n. (4.43) 

m^k 

As it is established in Sec 14. 3 1 the vector-function E(z) and the complex 
functions ipk( z ) express the same vector field in D k up to an isomorphism 
for E 0 = (—1,0). Therefore, application of the successive approximations to 
((4.43(1 yields the same series as (j3.3j) in D k - A detailed comparison can be 
performed. It shows, for instance, that the term Eo 
from (13.31) corresponds to (/ + pP kk ) *1 from (14.431) . 
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The solution of the problem (14.42j) depends on the shape of Dj~. Therefore, 
the first order approximation in v for the effective conductivity also depends 
on the shape and can be calculated by (14.3311 . If all the inclusions have the 
same form, the functions do not depend on k. Then, we arrive at the 

formula valid up to 0{y 2 ) 

A e ~ 1 + 2 pva, (4.44) 

where the shape factor a is introduced following [20], [34] 


a = 


Xi + ix 2 ) dxidx 2 . 


\DA 


'D i 


Application of the Pade (1, ^-approximation in v to ill. 1 II) yields 


, 1 + pva 

A e «--—. 

1 — pva 


(4.45) 


One can find limitations of the formulae il l. 1 Hi . i 1. 131) and a comparison of 
their accuracy in [34] . 

The higher order iterations in the cluster expansion are based on the 
integral equations on ^\z) 


k\ z ) + 2 L f Lk n l(t) ~z)dt = 

-pT.m^k 2 ^ifL m n2 m( t ) ~ z)dt + 1, 0 G D k (k = 1, 2, . . . , Tl) J 


1 , 2 ,- 


(4.46) 

The uniform convergence of the iterative scheme il 1. l(il) was proved in the case 
p = 1 (perfectly conducting inclusions) in Sec.4.9.2 of [35]. The higher order 
iterations (14 .46]) yield approximate analytical formulae for A e for arbitrary 
non-overlapping locations of different disks BJ.H0 and for ellipses 
pH]- An example of equal disks is discussed below in Sec 14.51 


4.5 Circular inclusions 

Let the inclusions be disks of radius r centered at the points a*.. Then, the 
normal unit vector nf.(t) to the circle \t — a^\ — r has the form rik(t) = t —^ ± . 
The inversion with respect to \t — a^l = r is defined by formula 

2 

z (k) = (4.47) 
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Consider the integral operator from 04.241) 


( Jkmf)(z ) 


1 

2ni 


7T7 / 


z)dt, z e D k , 


(4.48) 


where /(z) is analytic in the disk D m . This operator can be written in the 
form 


( Jkmf)(z) 



(4.49) 


where the function / is analytically continued into \z — a k \ > r. 

The function Ei(z) has poles of first order at the points z = miUi + 
m 2 C 02 (see [43]). The function ^ ._ f a j is analytic in the domain 

| z — a m | > r and the integral ( Jkmf){z ) can be calculated by residues at 
\z — a m \ > r. Introduce the operator 


M 

We have for m 


k 


--- I / f a m + = = ) 

z + miWi + m 2 uj 2 / \ * - a m + micui + m 2 uj 2 / 

(4.50) 


(Jkkf){z) 


E 


mi,m2EZ 


W m /) w 


(4.51) 


and for m ^ k 

<w)w = - E (wlWhw- < 4 - 52 ) 

mi,m2CZ 

The double series 04.50p ~ 04.52p are defined by the Eisenstein summation that 
correctly determines their convergence [6], [37] . 

Substitution of 04.50D - 04.52I) into 04.241) yields the system of functional 
equations 


N 


Vk{z) = 


m= 1 


mi ,7712 EZ 


(z) + 1, |z - a t \ < r (k = 1,2,.... JV), 


(4.53) 

where mi,m 2 run over integers in the sum )T)* with the excluded term 
m\ = m 2 = 0 for 77i = k. The functional equations 04.531) do not con¬ 
tain any integral. Each term ^1 Vml}m 2 '4 , m S j (z) can be treated as the shift 
operator written in the form of functional compositions. Such equations can 
be easily solved by use of symbolic computations. The difficulty related to 
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the infinite double summations in mi,m 2 can be overcome by application 
of the Eisenstein functions E n (z). Equations 04.531) can be solved by the 
contrast expansions in p (see Sec 14. 3 1) and by the cluster expansions in v (see 
Sec 14.41) . We refer to [36], [6], [7], [37] where this approach is explained in 
details. The computational efficiency of the method can be demonstrated by 
few analytical formulae below. Consider the cluster expansion in v = Nirr 2 
which is equivalent to the expansion in r 2 0 , m 

00 

^m(z) = ^2^\z)r 2n , \z-a k \ <r (k = 1,2, (4.54) 

n=0 

Few first functions in the expansion 04.54|) are given by the following exact 
formulae 

vl)m{z) = 1, ^m\z) = pY!u=i E ^ z ~ a fc)> 

'lpm(z) = P 2 Y!kM =1 E ^ ak - a ki) E 2 {z - a k ), 

(4.55) 

Tpm(z) = P 3 J2k,k u k 2 = 1 E ^ ak ~ a k 1 )E 2 (a kl - a k2 )E 2 (z - a k2 )~ 

2 P 2 Eu 1= i E s( a k ~ a kl )E 3 (z - a k ), 

where E n (z — a k ) is defined as E n (z — a k ) = E n (z — a k ) — (z — a k )~ p for k — m 
and E n ( 0) := E n ( 0). It is possible to proceed (14.551) and calculate the next 
approximations following the algorithm 0 , 0 . 

Let q be a natural number; k s runs over 1 to N, m q = 2,3,.... Let C 
denote the operator of complex conjugation. Introduce the multiple convo¬ 
lution sums 

e mi...m 5 — y [l+ipnH |_ mq )\ ^ 1 E mi { a ko — a k\ ) E m 2 ( a fci — a k 2 ) • • • {(l kq _ 1 — CL kq ). 

k 0 ki...k q 

(4.56) 

The integrals from (14.33[) are simplified by the mean value theorem for har¬ 
monic functions 

1 N 

An ~ i^i2 = 1 + 2 P u Jt y ^ J i’k(ak)- (4.57) 

V k =1 

Substitution of (14.551) and the next terms into (14.571) yields [31], [37] 

An — *Ai2 = 1 + 2 pi/(l + Aii/ + A 2 i/ 2 + ...), (4.58) 
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where 


— — 62 , A-2 — ^&22, A 3 — — 2/9 2 e33 + p 3 6222 ] j 

7T TT Z 7T 6 

= -7 [3p"e44 — 2 p 3 (e 332 + e 2 33 ) + p 4 e 2222 ] , 

7T* 

— ~7 [ — 4p 2 e55 + p 3 (3e442 + 6e 3 4 3 + 3e244) — 

7T° 

— 2p 4 (e 332 2 + e 233 2 + e 2 2 33 ) + p 5 e22222] , ( 4 . 59 ) 

— ~7 [ 5 p"e 66 — 4p 3 (e 25 5 + 3e 3 ,54 + 3e45 3 + 6552) + 

7T° 

+p 4 (3e2244 + 6e2 3 43 + 4e 3333 + 3e2442 + 6e 3 4 3 2 + 364422 ) — 

— 2p 5 (e22233 + ^22332 + ^23322 + e 33 222 ) + P & & 222222 )] • 

The next coefficients A n can be written in closed form by application of the 
algorithm presented in h, m- The deterministic values (14.56)) provide 
the benchmark to compute the effective conductivity of random composites. 
Following the Monte Carlo method one can take a representative sets of a e C 
to statistically calculate the expectations 

e mi ...m 9 (a)P(a)da. (4.60) 

The infinite set {e mi ... m ,rrij = 2,3,...} completely determine the random 
geometric structure of the considered class of composites and can be taken 
as the basic set in the RVE theory [30] . 

A fast algorithm to compute the convolution sums (14.56(1 is described in 
ra- It allows to compute A n (n < 20) within a given distribution of in 
few hours on an ordinary notebook when N = 64 and with the number of 
computational experiments 1500 in the Monte Carlo method. It follows from 
[9j that six terms (14.59)1 are sufficient to deduce an analytic formula for a 
uniform non-overlapping distribution of the disks. 

Formulae (14. 5811 - 04. 5911 corresponds to the cluster expansion. Application 
of the explicit iterations to the functional equations (14.53)) yields the contrast 
expansion as a series in p. One can see that the following third order formula 
in p takes place 

00 n —2 

An-iAi 2 = l+2pz/+2p 2 z/ 2 —+2p 3 z/ 3 V'(-l) n (n-l)e nn ^—+0(p 4 ). (4.61) 

7r z ' 7 T n 

n=2 

For macroscopically isotropic composites, the expectation of (14.61(1 has the 
form 

00 71 _2 

A e = 1 + 2 pu + 2pV + 2p 3 z/ 3 ^(-l) n (n - l)e nn ^—^- + 0(p 4 ), (4.62) 

71=2 



^-1 

a 4 

Ais 

^6 
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where the relation e 2 — r: from [33] is used. Another formula deduced in 
[33] can be useful in simulations and estimations of the third order term of 


gm>-gs2D 


(-ir 

jV n+1 


£ 

? 7 l=l 




(4.63) 


One can find the numerical values of e nn for uniform non-overlapping distri¬ 
butions in HD]. 

In order to compare the contrast expansion formula (14.62(1 and formula 
(I3.14P of the classic theory, we write (13.141) up to 0(p 4 ) 


pV(A e -l) 4 (A e + 1) = up - A^p 3 + 0(p 4 ). (4.64) 

Substitution of (14.621) into (14.641) yields 

( 00 n -2 \ 

^ ^ ~~ ^nn-^- ~ 1J • (4-65) 


The latter formula gives the Torquato-Milton parameter (see (20.59) and 
(20.66) from [42] ) 


Cl 



00 n —2 

^(-1) (n - l)e nn — 1 


(4.66) 


The three- and four- point contrast bounds on the effective conductivity 
are written by this parameter (see (21.33)-(21.35) and (21.42)-(21.44) from 

m)- 

Remark 4.4. Formula (I4.66|) is valid for arbitrary distributed circular inclu¬ 
sions. Any other shape of inclusions can be approximated by a disk packing. 

Remark 4.5. In accordance with (I4.37|) computation of the Torquato-Milton 
parameter Ci for general shapes of inclusions can be reduced to a sum of the 
triple integrals 



t 2 E l {t 2 - ti)Ei(ti - t)dtdtidt 2 


(4.67) 


Lk ^ -^ki ” ^- J k2 


that should be easier than computation of the corresponding triple integral 
(20.66) from [42] on the 3-point correlation function which is computed by 
another triple integral. The integral (I4.67|) can be calculated by residues for 
algebraic curves (see the case of ellipses discussed in |3l]). 
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5 Conclusion 


The iterative scheme (I4.37|) was constructively applied to special inclusions. 
An exact formula for regular array of disks (and solution of the problem 
(14.16(1 for N — 1) was obtained in [26] . [29] (see survey [37]). The term ”exact 
formula” should be explained precisely. Here, the exact formula means that 
A e is written as an expression which contains the geometrical parameters 
(the fundamental translation vectors ujj of the lattice Q, the radius of the 
disk) and the physical parameter p explicitly in symbolic form. The exact 
formula does not contains any parameter calculated by a numerical method 
(by an integral equation method or by truncation of an infinite system of 
equations). The exact formulae for A e in [37] was written in the form of the 
series similar to (13.10(1 in which all the coefficients C n were exactly written. 

Approximate analytical formulae for A e were obtained for arbitrary non¬ 
overlapping locations of different disks in [36] and of ellipses in [3T]. Here, the 
term ” approximate analytical formula” means that the effective conductivity 
A e was written in the form of the series 03.101) in which the first few coefficients 
C n were exactly written. Other results obtained by the GMS are discussed 
also in Sec l4.41 

In the previous works, the absolute convergence of the method was proved 
under geometrical restrictions [21]. Roughly speaking, these restrictions 
mean that each inclusion is far a way from others (see [42] and many pa¬ 
pers cited therein). However, it was proved in [35] that a modified method 
always uniformly converges, i.e., these geometrical restrictions are redundant. 
This is an interesting example of the difference between absolute and uniform 
convergence which shows that estimations on the absolute values or on the 
norm are too strong in comparison to the study of the uniform convergence 

[32]. 

The inclusions Dk can have exotic shapes that complicates the criterion of 
the dilute composite. This difficulty can be overcame by a conformal mapping 
of D onto a circular multiply connected domain. The capacity [14] (effective 
conductivity) is invariant under conformal mappings. Next, application of 
the functional equations for circular inclusions of different radii [36] can be 
applied to estimate their interactions. 

The GMS is effective in symbolic and numerical computations for non¬ 
overlapping inclusions of simple shapes (disks and ellipses) in the determin¬ 
istic statement. When an approximate analytical formula is deduced, the 
ensemble average is calculated by the straight forward computations. It is 
interesting to combine numerical methods effective for one-inclusion prob¬ 
lems with the analytically described interactions between inclusions by the 
GSM in Secl44l 
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Appendix. Convergence of the GMS 

This Appendix is devoted to study of the integral equations (I4.26p . (I4.29P 
in the space A. These equations can be shortly written in the form (I4.35p . 
We prove that equation (I4.35P has the unique solution for \p\ < 1 given 
by the series (14.361) obtained from (I4.35p by application of the successive 
approximations. 


29 














We use the following general result from [IS] - 

Theorem 6.1. Let B be a linear bounded operator in a Banach space B. If 
for any element f e B and for any complex number v satisfying the inequality 
\v\ < 1 equation 

'P = zaB'L + / (6.1) 

has a unique solution, then the unique solution of the equation 

'& = B'£ + f (6.2) 

can be found by the method of successive approximations. The approximations 
converge in B to the solution 


OO 

vf > = J2 Bk f■ (6-3) 

k =0 

Theorem 16.II can be applied to equation 'h = pAd> + 1 (see (14.35)1 ) in the 
Banach space A. 

Theorem 6.2. Let |p| < 1. Equation T = /n4T + 1 has a unique solu¬ 
tion. This solution can be found by the method of successive approximations 
convergent in the space A. 

Proof. Let \v\ < 1. Consider equations 

M Z ) = -VP El=l 2 h fLm n m (*) 'l>m(t)E 1 (t ~ z)dt + f k {z ), 2 G D k 

(k — 1,2,..., iV), 


Mt) = -VP [Whit) 1p k (t) + E m=1 2 h I’mitfEiit - ^)dtj + f k (t), 

t g L k (, k = 1,2,..., N), 
(6.4) 


where f k ^A k . 

Let 'ipk(z) be a solution of (14.29)1 . Introduce the function if(z) analytic in 
D and Holder continuous in its closure as follows 





<(i) i(^i(f 


z)dt, z G D. 


(6.5) 


One can see that if(z) is doubly periodic. The expression in the right hand 
part of (16.5)1 can be considered as Cauchy’s integral on torus 


$(*) 


1 

27 n 




z)dt 
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( 6 . 6 ) 













for z G Q\L with the density g,(t) = —upri^it) on L m . Using the 

property (I4.23P of (16. 6 p and (16.41) on L k we arrive at the M linear conjugation 
relation on each fixed curve L k 

kit) - f k {t) - -0(t) = -upn 2 k {t) il) k (t), t G L k . (6.7) 

Here, <f> + (t) = ipk{t) — f k (t ) and $~(f) = ip(t). Integration of (16.71) along L k 
yields 

= (j) k (t) - up <f> k (t) - g k {t) + 7 fc, t e L k (k = 1,2,...,1V), ( 6 . 8 ) 

where the functions 4>(z), (j) k {z) are analytic in .D, D k and belong to C^ 1,a ^ in 
the closures of the domains considered; g k (z) is a primitive function of f k (z ); 
7 ^ is an arbitrary constant. Moreover, the function <p(z) is quasi-periodic, 
hence 

4>{z) = <p(z) + cz, (6.9) 

where <p(z) is periodic and c is a constant. The constant 7 *, can be included 
into (pk{z) by the representation 7 k = ot k — up aj^. Then, (16. 8p becomes 

<p(t) = <p k (t) - up ip k (t) - g k (t) -ct, te L k , (k — 1,2,..., N), (6.10) 

It follows from |[8j and Appendix B of 123 that the general solution of the 
problem (16.101) has the following structure 

cp(z) = <p (1) (z) + (p (2) (z ), ip k (z) = <p { k \z) + <p£\z ), (6.11) 

where 

=<p ( k \t)-up(p ( k \t) - g k (t), t e L k , (k = 1,2,..., N) (6.12) 

and 

</? (2) (f) = y k \t) - up - ct, t e L k , (k = 1,2,.. .,N). (6.13) 

Each of the problems (16.121) and (16.131) (with fixed c) has a unique solution. 
The unique solution of the problem (16.131) can be easily found as ip' 2 \z) = 
cz and g^ k \z) = 0. Then, differentiation of (16. 1 11) (16. 1311 implies that the 
general solution of the problem (16.71) has the form 

ip(z) — 'ip^(z) + c, (6.14) 

where the function ijj k (z) is uniquely determined from the problem 

V> (1) 0) = il> k \t) + vpn 2 (t) - f k (t), t G L k (k = 1,2,..., N). (6.15) 
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(6.16) 


Comparison of (16.71) and (16.151) proves the theorem. 

The series 

OO 

T = s ^Tp k A k 1 

k =0 

converges absolutely for |p| < 1, since it can be represented in the form 
T = YlkLo s k PoA k l with sp 0 = p and |p| < |p 0 | < 1- The latter series con¬ 
verges absolutely with the rate |s| because the series YlT=o Po^ k ^- converges 
uniformly; hence, its general term tends to zero as k —> oo. 

Uniform convergence of (16.161) for \p\ = 1 can be justified by use of the 
arguments form [55] . 
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